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Abstract. Compact binary coalescence (CBC) is one of the most promising sources 
of gravitational waves. These sources are usually searched for with matched filters 
which require accurate calculation of the GW waveforms and generation of large 
template banks. We present a complementary search technique based on algorithms 
used in un-modeled searches. Initially designed for detection of un-modeled bursts, 
which can span a very large set of waveform morphologies, the search algorithm 
presented here is constrained for targeted detection of the smaller subset of CBC 
signals. The constraint is based on the assumption of elliptical polarisation for signals 
received at the detector. We expect that the algorithm is sensitive to CBC signals in a 
wide range of masses, mass ratios, and spin parameters. In preparation for the analysis 
of data from the fifth LIGO- Virgo science run (S5), we performed preliminary studies 
of the algorithm on test data. We present the sensitivity of the search to different 
types of simulated binary black hole waveforms. Also, we discuss how to extend the 
results of the test run into a search over all of the current LIGO- Virgo data set. 
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1. Introduction 

The merging of massive compact binary systems is one of the most promising sources 
for detection of gravitational waves with interferometric detectors. Typically, these 
sources are composed of binary systems of neutron stars and/or black holes, and the 
gravitational waves radiated from these sources are coupled to new physics. These 
sources are critical for a test of general relativity, and provide a valuable way to study 
the strong-field regime. Their orbits slowly decay because energy is lost in the emission of 
gravitational radiation, and the frequency of the emitted quasi-sinusoidal gravitational 
waves increases over time. This semi-periodic stage is called the inspiral. Due to the 
decay of the orbit, the pair eventually comes close enough for strong-field effects to 
dominate and no stable orbit is available to the system. The two compact objects 
then merge together to form a single black hole, releasing a burst of gravitational wave 
radiation. This stage is referred to as the merger stage. Finally, in the ringdown 
stage, the newly formed single black hole spins down and radiates gravitational waves 
in damped quasi-normal modes. The more massive a binary system is, the more energy 
is expected to be radiated in gravitational waves. Higher mass also corresponds to a 
lower frequency bandwidth of signal content. These frequencies fall within the most 
sensitive region for ground based gravitational wave interferometers for stellar mass and 
some intermediate mass compact binaries. 

A number of searches for gravitational waves from compact binary sources have 
been undertaken in the past, some of which are similar to the search presented 
here[T]. However, most searches[2], [3l H] are based on the matched filter technique 
which uses templates for detection of signals of known waveforms [5]. These templates 
span a parameter space: a set of source parameters such as total mass, mass ratio, 
spin, and other properties which describe the binary system. It is known that a 
matched-filter approach is optimal, however, this approach requires the construction 
of large template banks using accurate theoretical models which is not always possible. 
Currently, a template bank which covers all possible gravitational wave signals from 
binary coalescences does not exist. For example, it is expected that some binary mergers 
will involve components which are spinning[6]. It has been suggested [TJ [8] that coalescing 
binaries with spinning components will have strong effects on the merger stage of the 
waveform for masses greater than 50 MQbecause of the contribution of higher order 
modes. While some progress has been made in constructing template banks for the 
inspiral stage of spinning binaries [9| [TO], there are no waveforms describing all stages of 
coalescence as well as being currently available for use in template based searches. 

Despite problems in analytically formulating the late stages of binary coalescences, 
recent progress in numerical relativity has allowed for a calculation of the waves emitted 
during the merger stage of two binary black holes. As a result, template banks have 
been supplemented by numerical simulations of the full inspiral and merger stages of 
the black hole binary system. This has enabled the construction of high mass (> 35Mq) 
template banks used in searches for those binaries. There are still difficulties involved 
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in searching for high mass binaries. The calculation of the numerical waveforms is 
computationally costly. The generation process makes use of numerical waveforms, but, 
because of the cost, the waveforms produced are still only approximate in the merger 
region. 

In addition to template based searches, there is another, more general, class of 
searches that have been carried out for gravitational waves from un-modeled, short 
duration sources. They are sensitive to a wide class of waveform morphologies and do 
not require template banks. The detection of un-modeled, or "burst" , sources is based 
on excess signal power and correlations of signals detected in different gravitational wave 
detectors. Because an un-modeled search does not require templates, it is useful in the 
regions of the compact binary parameter space where there are no waveforms available 
for use or where the waveforms are difficult to compute. 

In this paper we present a robust method for detection of gravitational waves from 
compact binary coalescences. This method is expected to be less sensitive than the 
optimal matched filters, but encompasses a wider class of compact binary sources that 
may be missed by the template searches. 

2. Search Method 

We first present a brief description of the Coherent WaveBurst algorithm used for 
searches of LIGO data for gravitational waves from sources which may not be modeled. 
For a detailed discussion we refer to its implementation details in the references 
After the discussion of its un-modeled stage of operation, we will describe a new 
constraint introduced into the search in the context of detection of binary coalescence 
signals. We will also outline the execution of a test run using a week of simulated 
noise data and present the results. These results will include the injection of simulated 
waveforms into the aforementioned data set, and use those results to estimate the 
sensitivity of the search to compact binary sources. 

2.1. Coherent WaveBurst 

Coherent WaveBurst is a coherent network algorithm which searches for un-modeled, 
short duration bursts of gravitational waves in detector strain data using information 
from all detectors in a network simultaneously [T2] . Therefore, events are reconstructed 
for the network, not individual detectors. This is compared to generating detector 
specific events and attempting to coincide them across detectors a posteriori^^. The 
algorithm has been used in the fourth [T^ and fifth[15] LIGO science runs in searches for 
gravitational waves from potential sources such as supernovas, gamma ray bursts, and 
others. Because of the lack of source models, a statistical approach is employed which 
is based mainly on excess cross-correlated power inside the detector network [T7] . 

In practice, the method described above is performed using time-frequency 
decomposition. Detector data is read into Coherent WaveBurst and decomposed to 
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multiple time-frequency resolutions. This decomposition transforms a sampled time 
series into a time-frequency map, with the values of locations on the map representing 
the energy amplitude in that time-frequency location. Each detector has its own time- 
frequency map, and in order to create coherent events, the detectors' time series are 
delayed relative to each other with respect to a set of possible travel times of the 
gravitational wave. The maps are then combined together for each delay. 

Coherent WaveBurst identifies event candidates as regions on the time-frequency 
plane which have excess power. A clustering procedure is employed to identify the same 
time-frequency event across multiple decomposition resolutions. Once all of the clusters 
are identified, then coherent statistics are calculated for each event candidate. 

2.2. Un-modeled Likelihood 

Coherent WaveBurst conducts the above method in the absence of a source model by use 
of the likelihood ratio [H], which has been studied before [T8l [19]. A similar framework 
for detection of the inspiral stage waveforms using the likelihood approach has also 
been suggested [20]. but in the context of a template based search. We describe here the 
construction of a likelihood ratio functional in the model independent case assuming 
stationary and gaussian noise. We consider a generic waveform h composed of two 
polarisation states h+ and hx ■ The effect of a gravitational wave on an interferometer 
is characterized by a detector response ^ 



The detector response is a function of the sky coordinates {9, (p) through the antenna 
patterns (F+, F^), which represent the interferometer's sensitivity to the respective 
polarisation states. In Coherent WaveBurst, the likelihood functional is the ratio of the 
joint probability P{x\h) of a gravitational wave signal being present in the data to the 
joint probability P{x\0) of no gravitational wave signal. The ratio is given by: 



In the above equation Xkj represents the sampled output strain from an individual 
detector k and time- frequency location j. The denominator, alj, is the noise variance in 
the kth detector. This is summed over each of the K detectors in the network and time- 
frequency region composed of time-frequency samples. This functional is analytically 
maximized to determine the estimators of /i+ and amplitudes for each sample of the 
data. We consider these amplitudes to be free parameters used in the variation of the 
likelihood functional. With two (possibly independent) polarisation states, this allows 
at most 2N degrees of freedom to fit with K x N samples from the data. Because the 
detector responses are also a function of the sky coordinates via the antenna patterns, 
the likelihood functional is also maximized over a set of sky locations. The location 
at the maximum is an estimate of the source location in the sky. In practice, this 
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process involves a large number of free parameters to maximize simultaneously, and the 
likelihood functional may allow for unphysical solutions for the reconstructed waveforms. 
In order to mitigate these problems, constraints on the likelihood functional are used 
to suppress unwanted solutions. In the search presented, we are also using constraints 
to introduce the model dependence as described in section 12.31 The intended purpose 
of this constraint is to increase sensitivity to compact binary sources as well as better 
rejection of noise induced events in real detector data. 

2.3. Elliptical Constraint 

Since we intend to search for compact binary signals with more definite physical 
properties than in the un-modeled case, it is possible to further constrain the likelihood 
functional to better reflect prior knowledge about the signal model. As an inspiralling 
compact object is a rotating system, the waves produced should be elliptically polarised, 
and therefore our goal is to enhance the reconstruction of elliptically polarised signals. 
The enforcement of this constraint should also cause incoherently polarised signals to 
be poorly reconstructed, which is important in the rejection of glitches. This constraint 
is expressed by modifying the meaning of the estimators for /i_|_and h^: 

h+ ^ h 
hy. ^ h 

= F+{ij)h + F,{ij)ah. (3) 

By this modification, ^' now represents the detector response to an elliptically 
polarised gravitational wave composed of h and its ninety degrees phase shifted 
counterpart h. We will denote all quantities with this phase shift with a tilde. The 
principal effect of the additional constraint on the algorithm is to reduce number of 
free parameters required for waveform reconstruction. In the model independent case, 
there are two independent polarisation states of the gravitational wave. For a set of N 
sample measurements in the time series, this corresponds to 2N free parameters that 
the algorithm must fit. However, elliptical polarisation demands that one polarisation 
state h is related to the other /i by a phase shift of ninety degrees. We can exploit this 
relationship because effectively only one polarisation state is required to be fit, as its 
phase shifted polarisation state should be identical up to an amplitude factor. The use 
of this relationship reduces the number of free parameters required to for h and two 
additional parameters i/j and a. The plane that the states h and h lie in are rotated with 
an angle ifj with respect to the detector plane. This angle is the polarisation angle of the 
wave. Also, the ellipticity parameter a is introduced. This parameter is proportional to 
the relative strength of the polarisation states where the unity cases (a = ±1) correspond 
to circular polarisation and the zero case to linear polarisation. For binary sources, this 
parameter is implicitly related to the inclination angle of the source with respect to the 
observer. 
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Coherent WaveBurst makes use of the constraint by substituting ^' in the definition 
of the hkehhood ratio in ([2]). The estimators for h and h are then obtained in a similar 
varying procedure as in the un-modeled case. 

2.4- Selection Criteria 

Usually data collected by gravitational wave detectors is non-stationary and polluted 
with environmental and instrumental artifacts. Searches use selection criteria to 
distinguish between gravitational wave signals and "glitches" in the data from those 
artifacts. In previous searches [H] using Coherent WaveBurst, there were a set of four 
statistics, calculated within the likelihood approach, which were chosen for selection of 
coherent events. We presently describe those statistics, but refer to [H] for a fuller 
description. 

We define the network correlation coefficient cc, which tests the overall consistency 
of the event candidate based on the comparison of reconstructed correlated signal energy 
Ec and residual noise energy En. The correlated signal energy is a weighted sum of 
the correlated energies from all pairs of detectors in the network, and the residual 
noise energy represents the energy that remains after the reconstructed signal energy is 
subtracted from the total energy of the event. 



This criteria is intended to exploit the poor reconstruction of glitches. A glitch often 
can have large residual noise energy as compared to its correlated energy. The network 
correlation coefficient varies between zero and one, where true gravitational waves should 
have a value of cc near unity, and glitches will have values of cc significantly less than 
unity. 

We define a penalty factor Pj, which is a measure of consistency of reconstructed 
detector responses. This is a selection criteria for event candidates which have 
unphysically reconstructed responses. Such candidates have a reconstructed signal 
energy which is greater than the total energy in a given detector. This statistic is 
defined by: 



The quantities {x\) and (^|) are the summed data sample amplitudes and 
reconstructed detector responses for a selected time-frequency region in an individual 
detector k. If an event is inconsistent, the detector which is maximally inconsistent is 
chosen. If the event is consistent {{^1) < {x\) for all A;), then this parameter is defined 
to be unity. 

The next criteria is the energy disbalance A, which is similar to Pf but indicates the 
overall mismatch between the reconstructed amplitudes and the amplitudes recorded by 



Pf = min 

k 




(5) 



Burst CBC search 



Table 1: Overview of the parameters used in the injection simulations. 





EOBNR.vla 


EOBNR.vlb 


total mass (Mq) 


5 < Mtotal < 100 


100 < Mt^tai < 500 


component masses (Mq) 


1 < M, < 95 


40 < M, < 400 


minimum mass ratio 


0.012 


0.1 



the detectors: 



A 



(6) 



It characterizes by the mismatch between the reconstructed energy of the event and 
the energy of the data. A glitch in one detector can lead to large detector responses 
being reconstructed in other detectors. For example, the (^f ) term will exceed the {xk^k) 
term and produce large values of A. Ideally, this statistic should be zero. 

All event candidates are assigned a coherent network amplitude. This statistic is 
the main detection statistic and used to rank events: 



V 



lEcCC 



(7) 



Because glitches are characterized by a small correlated energy, and therefore can 
also have a small network correlation coefficient, this statistic is used together with cc 
to reject various types of glitches. The combination of the two statistics allows the 
rejection of both consistent but weak glitches, and strong but less consistent glitches. 
It is used as a detection statistic for event candidates after all other thresholds have 
been applied. The threshold of rj used is a function of false alarm rate selected for that 
analysis. 



3. Simulation and Estimation of the Detection Efficiency 

In order to estimate the sensitivity of the search to gravitational waves from compact 
binary sources, we make use of the EOBNR (Effective One Body Numerical Relativity) 
family of waveforms [2 ![ [22| [23] . In this model, analytical inspiral- plunge waveforms 
computed in the Effective One Body formalism |24j are calibrated to the numerical 
relativity simulations, and matched to the quasi-normal modes of the ringdown to 
produce complete coalescence signals. The waveforms are only constructed for case 
of non-spinning black hole binaries. 

For estimation of the detection efficiency these waveforms are used by injecting 
them into the detector strain data, and then the algorithm is applied to detect and 
reconstruct them. The injections cover a wide range of source parameters. Table [1] 
describes the two sets of injections tested, delineated by total mass. Common between 
the injection sets were uniform distributions in sky position, inclination angle, and time 



Burst CBC search 



8 



Table 2: Values of selection criteria used in the analysis. 



network correlation (cc) 


> 


0.6 


penalty (P/) 


> 


0.6 


energy disbalance (A) 


< 


0.35 


coherent network ampltide (77) 


> 


3.0 



between injection. All of the waveforms were generated at a distance of 1 Mpc and 
scaled to twenty other distances between 100 kpc and 1 Gpc. 

The simulated noise data into which simulated signals are injected has a colored 
gaussian profile which approximates the LIGO network during its fifth science run. A set 
of background events are generated from the simulated detector data with no simulated 
signal injections. These background events are created by shifting detectors in time with 
respect to each other beyond the maximum possible light travel time of a gravitational 
wave through Earth. This background set is used to tune thresholds and estimate 
the false detection rate due to random coincidences by identifying loud background 
events and fixing the selection criteria thresholds as found in table [2] to exclude these 
events from the analysis. The penalty factor, network correlation coefficient, and energy 
disbalance thresholds are applied to reject loud but inconsistent background events, and 
the false detection rate (10^^ Hz for this analysis) picked then determines the threshold 
on weak background events via the coherent network amplitude r]. 

Once this stage is complete, then the simulated compact binary waveforms are 
injected into the data, and the algorithm is applied to detect these. Simulated event 
candidates are passed through the same set of selection criteria which were applied to 
the background. 

3.1. Simulation Detection Results 

Estimated detection efficiencies are shown for both sets of waveforms in figure [H From 
figure [ic] and figure [Tdl it is shown that we are sensitive to the wide range of parameters 
exhibited by these injection sets. In particular, our detection efficiencies at close 
distances approach unity, and therefore we are sensitive to the entire set of mass ratios 
used in the injection sets. It is also evident from figure [Ia| and figure [lb] that we are 
sensitive to a wide range of total system mass, enabling investigation of both stellar 
mass black hole binaries and some intermediate mass binary black hole systems. Also, 
almost all of the total mass bins in figures [Tal and [Tbl have a 90% detection efficiency up 
to 20 Mpc. This distance encloses the Virgo cluster, an abundant collection of potential 
hosts for sources. 



Burst CBC search 



9 




Distance (Mpc) Distance (Mpc) 

(c) (d) 

Figure 1: Network detection efficiencies for the low mass (left) and high mass (right) 
sets. The top plot shows efficiency over the total mass and injection distance of the 
system while the bottom plot is the efficiency after integration over the total mass 
dimension. 

3.2. Estimation of Visible Luminosity 

We estimate the number of hosts visible to the search, via their luminosity, by use 
of the network detection efficiencies and a catalog of galaxies. The catalog compiled 
by the compact binary coalescence group in LIGO[25] with luminosity entries is used 
because of the unknown relationship between compact binary source populations and 
their host galaxies. A simple relationship between the stellar birth rate, as measured 
by the blue light luminosity of a galaxy, and the binary coalescence rate is assumed. 
Because of this simplification, luminosities calculated by compact binary searches are 
typically calculated in units of blue light luminosity Liq. This luminosity represents the 
output of 10^° G type stars in the blue light (445 nm) optical ffiter[26]. For reference, the 
Milky Way galaxy is approximately 1.7 Liq. The total luminosity visible to a search, as 
measured in units Lio, is called the cumulative luminosity. The cumulative luminosity 
is a measure of the effective number of galaxies surveyed by a search and is used in rate 
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calculations in template based searches. We calculate the efficiency e of the network 
and use this to convolve with each galaxy having a location {ri,6i,(f)i) and luminosity 
measurement Lf. 

CL = Y.<^^^^^^^^)L^■ (8) 

i 

In this study, we have considered a wider range of masses than previously examined 
in other searches, and as a consequence the increased gravitational wave energy output 
for higher mass binaries and the corresponding increase in search volume available allows 
for an order of magnitude increase in detected luminosity. In practice, the cost of a 
widening the mass range is reflected in the subsequent requirement of less restrictive 
selection criteria, which leads to a higher false detection rate. However, the simulated 
nature of the data does not contain those transients, which tend to form the outliers in 
the distribution of 77, and we are forced to disregard the effects on the false detection 
rate. As a consequence of using simulated data, the threshold on 77 is decreased and 
will therefore lead to comparatively better efficiencies. Matched ffiter searches using the 
most current set of LIGO data obtained a cumulative luminosity of 5500 Lio for binary 
black hole systems of total mass less than 35 M0 14] . In our test with simulated data, we 
estimate a cumulative luminosity of above 22000 Lio for the low mass injection set. The 
matched filter cumulative luminosities quoted are for searches involving the Livingston 
detector and both Hanford detectors. 

4. Future Work 

There are a number of items that have not been addressed in the test run presented, 
but will be included in the full search. Because the test run used simulated noise, rather 
than actual detector data, we expect that some of these estimates are optimistic. There 
are a number of issues involved in actual detector data. One of the more difficult issues 
to deal with are instrumental and environmental transients falsely identified as event 
candidates, however, the enhancement of our selection procedures with regards to the 
elliptical constraint should be able to mitigate this problem. More study should allow 
us to characterize how well we can reject glitches in real detector data, and therefore 
how effective the elliptical constraint is in the reduction of glitches. 

The latest data taking run by LIGO has produced the equivalent of a year of 
data in which all three interferometers were taking science quality data simultaneously. 
Additionally, in the later part of this run, the Virgo interferometer was also in 
coincidence with the LIGO network and the data taken from the Virgo run will 
be analyzed jointly with LIGO data. It is also possible to perform more extensive 
simulations with these and other waveforms |2 7] in order to place plausible upper limits 
on the rate of compact binary coalescences in the universe in the event that gravitational 
wave signals are not confidently detected. 

It is expected that LIGO and Virgo will start a joint run with upgraded detectors 
before the end of 2009. The search presented can also be used for the analysis of this 
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new, more sensitive data. 
5. Summary 

We have presented some of the motivations for carrying out a search for compact binary 
coalescences with an un-modeled technique. These motivations include detecting events 
which may not be well modeled analytically and therefore may not be detected by 
a template based search. To this end, we have outlined the plans for a search for 
binary coalescences using a weakly-modeled search algorithm which has been modified 
to better detect the intended source model of the search. In this regard, the search is 
intended to be a complement to existing template based searches. These modifications 
are intended to improve the reconstruction of elliptically polarised gravitational wave 
signals. We have also presented the results of a test run using simulated LIGO data and 
simulated binary coalescence waveforms in order to estimate the efficiency of the search 
in detection of binary coalescence signals. We have demonstrated the ability to detect 
90% of injected signals at 10 Mpc for systems of total mass between 100-500 MQand 20 
Mpc for systems of 25-100 Mq. 
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